Deletion in KRT75L4 linked to frizzle feather in Xiushui Yellow Chickens.

Bird feathers are the product of interactions between natural and artificial selection. Feather-related traits are important for chicken selection and breeding. Frizzle feather is characterized by the abnormally development of feathers in chickens. In the current study, frizzle feather characteristics were observed in a local breed called Xiushui Yellow Chicken in Jiangxi, China. To determine the molecular mechanisms that underlie frizzle feather in Xiushui Yellow Chicken, four populations of three breeds (Xiushui Yellow Chicken with frizzle feathers, Xiushui Yellow Chicken with normal feathers, Guangfeng White-Ear Yellow Chicken, and Ningdu Yellow Chicken) were selected for whole-genome resequencing. Using a comparative genome strategy and genome-wide association study, a missense mutation (g.5281494A>G) and a 15-bp deletion (g.5285437-5285451delGATGCCGGCAGGACG) in KRT75L4 were identified as candidate mutations associated with frizzle feather in Xiushui Yellow Chicken. Based on genotyping performed in a large Xiushui Yellow Chicken population, the g.5285437-5285451delGATGCCGGCAGGACG mutation in KRT75L4 was confirmed as the putative causative mutation of frizzle feather. These results deepen the understanding of the molecular mechanisms responsible for frizzle feather, as well as facilitating the molecular detection and selection of the feather phenotype in Xiushui Yellow Chickens.

OTOGL, a gelforming mucin protein, is nonessential for male germ cell development and spermatogenesis in mice.

Otogelin-like protein (encoded by Otogl) was highly structural similar to the gelforming mucin proteins. Although human OTOG mutations have been linked to deafness, the biological function of OTOGL in male germ cell development remains enigmatic. In screening 336 patients with non-obstructive azoospermia (NOA), OTOGL displays the high mutant ratio (13.99 %). Then, we examined the expression of OTOGL in developing mouse testes. Otogl mRNA and protein are continually expressed in postnatal developing testes from postnatal day 0 (P0) testes to P21 testes exhibiting a decreased trend with the age growth. We thus generated a global Otogl knockout mouse (KO) model using the CRISPR/Cas9 technology; however, Otogl KO mice displayed normal development and fertility. Further histological analysis of Otogl knockout mouse testes revealed that all types of spermatogenic cells are present in Otogl KO seminiferous tubules. Together, our study suggested that OTOGL is nonessential for male germ cell development and spermatogenesis.

Rescue of two trafficking-defective variants of the neuronal glycine transporter GlyT2 associated to hyperekplexia.

Hyperekplexia is a rare sensorimotor syndrome characterized by pathological startle reflex in response to unexpected trivial stimuli for which there is no specific treatment. Neonates suffer from hypertonia and are at high risk of sudden death due to apnea episodes. Mutations in the human SLC6A5 gene encoding the neuronal glycine transporter GlyT2 may disrupt the inhibitory glycinergic neurotransmission and cause a presynaptic form of the disease. The phenotype of missense mutations giving rise to protein misfolding but maintaining residual activity could be rescued by facilitating folding or intracellular trafficking. In this report, we characterized the trafficking properties of two mutants associated with hyperekplexia (A277T and Y707C, rat numbering). Transporter molecules were partially retained in the endoplasmic reticulum showing increased interaction with the endoplasmic reticulum chaperone calnexin. One transporter variant had export difficulties and increased ubiquitination levels, suggestive of enhanced endoplasmic reticulum-associated degradation. However, the two mutant transporters were amenable to correction by calnexin overexpression. Within the search for compounds capable of rescuing mutant phenotypes, we found that the arachidonic acid derivative N-arachidonoyl glycine can rescue the trafficking defects of the two variants in heterologous cells and rat brain cortical neurons. N-arachidonoyl glycine improves the endoplasmic reticulum output by reducing the interaction transporter/calnexin, increasing membrane expression and improving transport activity in a comparable way as the well-established chemical chaperone 4-phenyl-butyrate. This work identifies N-arachidonoyl glycine as a promising compound with potential for hyperekplexia therapy.

Pathogenic missense protein variants affect different functional pathways and proteomic features than healthy population variants.

Missense variants are present amongst the healthy population, but some of them are causative of human diseases. A classification of variants associated with "healthy" or "diseased" states is therefore not always straightforward. A deeper understanding of the nature of missense variants in health and disease, the cellular processes they may affect, and the general molecular principles which underlie these differences is essential to offer mechanistic explanations of the true impact of pathogenic variants. Here, we have formalised a statistical framework which enables robust probabilistic quantification of variant enrichment across full-length proteins, their domains, and 3D structure-defined regions. Using this framework, we validate and extend previously reported trends of variant enrichment in different protein structural regions (surface/core/interface). By examining the association of variant enrichment with available functional pathways and transcriptomic and proteomic (protein half-life, thermal stability, abundance) data, we have mined a rich set of molecular features which distinguish between pathogenic and population variants: Pathogenic variants mainly affect proteins involved in cell proliferation and nucleotide processing and are enriched in more abundant proteins. Additionally, rare population variants display features closer to common than pathogenic variants. We validate the association between these molecular features and variant pathogenicity by comparing against existing in silico variant impact annotations. This study provides molecular details into how different proteins exhibit resilience and/or sensitivity towards missense variants and provides the rationale to prioritise variant-enriched proteins and protein domains for therapeutic targeting and development. The ZoomVar database, which we created for this study, is available at fraternalilab.kcl.ac.uk/ZoomVar. It allows users to programmatically annotate missense variants with protein structural information and to calculate variant enrichment in different protein structural regions.

Barrier-to-autointegration factor: a first responder for repair of nuclear ruptures.

The nuclear envelope (NE) is a critical barrier between the cytosol and nucleus that is key for compartmentalization within the cell and serves an essential role in organizing and protecting genomic DNA. Rupturing of the NE through loss of constitutive NE proteins and/or mechanical force applied to the nucleus results in the unregulated mixing of cytosolic and nuclear compartments, leading to DNA damage and genomic instability. Nuclear rupture has recently gained interest as a mechanism that may participate in various NE-associated diseases as well as cancer. Remarkably, these rupturing events are often transient, with cells being capable of rapidly repairing nuclear ruptures. Recently, we identified Barrier-to-Autointegration Factor (BAF), a DNA-binding protein involved in post-mitotic NE reformation and cytosolic viral regulation, as an essential protein for nuclear rupture repair. During interphase, the highly mobile cytosolic BAF is primed to monitor for a compromised NE by rapidly binding to newly exposed nuclear DNA and subsequently recruiting the factors necessary for NE repair. This review highlights the recent findings of BAF's roles in rupture repair, and offers perspectives on how regulatory factors that control BAF activity may potentially alter the cellular response to nuclear ruptures and how BAF may participate in human disease.

Missense mutation of Fmr1 results in impaired AMPAR-mediated plasticity and socio-cognitive deficits in mice.

Fragile X syndrome (FXS) is the most frequent form of inherited intellectual disability and the best-described monogenic cause of autism. CGG-repeat expansion in the FMR1 gene leads to FMR1 silencing, loss-of-expression of the Fragile X Mental Retardation Protein (FMRP), and is a common cause of FXS. Missense mutations in the FMR1 gene were also identified in FXS patients, including the recurrent FMRP-R138Q mutation. To investigate the mechanisms underlying FXS caused by this mutation, we generated a knock-in mouse model (Fmr1R138Q) expressing the FMRP-R138Q protein. We demonstrate that, in the hippocampus of the Fmr1R138Q mice, neurons show an increased spine density associated with synaptic ultrastructural defects and increased AMPA receptor-surface expression. Combining biochemical assays, high-resolution imaging, electrophysiological recordings, and behavioural testing, we also show that the R138Q mutation results in impaired hippocampal long-term potentiation and socio-cognitive deficits in mice. These findings reveal the functional impact of the FMRP-R138Q mutation in a mouse model of FXS.

A Missense Mutation in a Large Subunit of Ribonucleotide Reductase Confers Temperature-Gated Tassel Formation.

Temperature is a major factor regulating plant growth. To reproduce at extreme temperatures, plants must develop normal reproductive organs when exposed to temperature changes. However, little is known about the underlying molecular mechanisms. Here, we identified the maize (Zea mays) mutant thermosensitive vanishing tassel1-R (tvt1-R), which lacks tassels at high (restrictive) temperatures due to shoot apical meristem (SAM) arrest, but forms normal tassels at moderate (permissive) temperatures. The critical stage for phenotypic conversion in tvt1-R mutants is V2 to V6 (Vn, where "n" is the number of leaves with collars visible). Positional cloning and allelism and complementation tests revealed that a G-to-A mutation causing a Arg277-to-His277 substitution in ZmRNRL1, a ribonucleotide reductase (RNR) large subunit (RNRL), confers the tvt1-R mutant phenotype. RNR regulates the rate of deoxyribonucleoside triphosphate (dNTP) production for DNA replication and damage repair. By expression, yeast two-hybrid, RNA sequencing, and flow cytometric analyses, we found that ZmRNRL1-tvt1-R failed to interact with all three RNR small subunits at 34°C due to the Arg277-to-His277 substitution, which could impede RNR holoenzyme (α2ß2) formation, thereby decreasing the dNTP supply for DNA replication. Decreased dNTP supply may be especially severe for the SAM that requires a continuous, sufficient dNTP supply for rapid division, as demonstrated by the SAM arrest and tassel absence in tvt1-R mutants at restrictive temperatures. Our study reveals a novel mechanism of temperature-gated tassel formation in maize and provides insight into the role of RNRL in SAM maintenance.

Comprehensive characterization of amino acid positions in protein structures reveals molecular effect of missense variants.

Interpretation of the colossal number of genetic variants identified from sequencing applications is one of the major bottlenecks in clinical genetics, with the inference of the effect of amino acid-substituting missense variations on protein structure and function being especially challenging. Here we characterize the three-dimensional (3D) amino acid positions affected in pathogenic and population variants from 1,330 disease-associated genes using over 14,000 experimentally solved human protein structures. By measuring the statistical burden of variations (i.e., point mutations) from all genes on 40 3D protein features, accounting for the structural, chemical, and functional context of the variations' positions, we identify features that are generally associated with pathogenic and population missense variants. We then perform the same amino acid-level analysis individually for 24 protein functional classes, which reveals unique characteristics of the positions of the altered amino acids: We observe up to 46% divergence of the class-specific features from the general characteristics obtained by the analysis on all genes, which is consistent with the structural diversity of essential regions across different protein classes. We demonstrate that the function-specific 3D features of the variants match the readouts of mutagenesis experiments for BRCA1 and PTEN, and positively correlate with an independent set of clinically interpreted pathogenic and benign missense variants. Finally, we make our results available through a web server to foster accessibility and downstream research. Our findings represent a crucial step toward translational genetics, from highlighting the impact of mutations on protein structure to rationalizing the variants' pathogenicity in terms of the perturbed molecular mechanisms.

Studying the Effects of ACE2 Mutations on the Stability, Dynamics, and Dissociation Process of SARS-CoV-2 S1/hACE2 Complexes.

A highly infectious coronavirus, SARS-CoV-2, has spread in many countries. This virus recognizes its receptor, angiotensin-converting enzyme 2 (ACE2), using the receptor binding domain of its spike protein subunit S1. Many missense mutations are reported in various human populations for the ACE2 gene. In the current study, we predict the affinity of many ACE2 variants for binding to S1 protein using different computational approaches. The dissociation process of S1 from some variants of ACE2 is studied in the current work by molecular dynamics approaches. We study the relation between structural dynamics of ACE2 in closed and open states and its affinity for S1 protein of SARS-CoV-2.

Collateral fitness effects of mutations.

The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 ß-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.

An Epilepsy-Associated SV2A Mutation Disrupts Synaptotagmin-1 Expression and Activity-Dependent Trafficking.

The epilepsy-linked gene SV2A, has a number of potential roles in the synaptic vesicle (SV) life cycle. However, how loss of SV2A function translates into presynaptic dysfunction and ultimately seizure activity is still undetermined. In this study, we examined whether the first SV2A mutation identified in human disease (R383Q) could provide information regarding which SV2A-dependent events are critical in the translation to epilepsy. We utilized a molecular replacement strategy in which exogenous SV2A was expressed in mouse neuronal cultures of either sex, which had been depleted of endogenous SV2A to mimic the homozygous human condition. We found that the R383Q mutation resulted in a mislocalization of SV2A from SVs to the plasma membrane, but had no effect on its activity-dependent trafficking. This SV2A mutant displayed reduced mobility when stranded on the plasma membrane and reduced binding to its interaction partner synaptotagmin-1 (Syt1). Furthermore, the R383Q mutant failed to rescue reduced expression and dysfunctional activity-dependent trafficking of Syt1 in the absence of endogenous SV2A. This suggests that the inability to control Syt1 expression and trafficking at the presynapse may be key in the transition from loss of SV2A function to seizure activity.SIGNIFICANCE STATEMENT SV2A is a synaptic vesicle (SV) protein, the absence or dysfunction of which is linked to epilepsy. However, the series of molecular events that result in this neurological disorder is still undetermined. We demonstrate here that the first human mutation in SV2A identified in an individual with epilepsy displays reduced binding to synaptotagmin-1 (Syt1), an SV protein essential for synchronous neurotransmitter release. Furthermore, this mutant cannot correct alterations in both Syt1 expression and trafficking when expressed in the absence of endogenous SV2A (to mimic the homozygous human condition). This suggests that the inability to control Syt1 expression and trafficking may be key in the transition from loss of SV2A function to seizure activity.

Importance of asparagine-381 and arginine-487 for substrate recognition in CYP4Z1.

The human cytochrome P450 enzyme CYP4Z1 remains an understudied enzyme despite its association with poor prognosis and overexpression in breast cancer. Hence, CYP4Z1 has previously been suggested as an anti-breast cancer target. In the present study we employed extended mutation analysis to increase our understanding of the substrate binding mode of this enzyme. In a combined in vitro and in silico approach we show for the first time that residue Arg487 plays an important role in substrate recognition and binding of CYP4Z1. Using a large array of recombinant CYP4Z1 mutants we show that, apart from Asn381, all other postulated binding residues only play an auxiliary role in substrate recognition and binding. Different substrate interaction motifs were identified via dynamic pharmacophores (dynophores) and their impact on catalytically competent substrate binding was classified. These new insights on the substrate recognition and binding mode represent an important step towards the rational design of CYP4Z1 prodrugs and guide further investigations into the so far poorly understood physiological role of CYP4Z1.

The origins of dengue and chikungunya viruses in Ecuador following increased migration from Venezuela and Colombia.

BACKGROUND: In recent years, Ecuador and other South American countries have experienced an increase in arboviral diseases. A rise in dengue infections was followed by introductions of chikungunya and Zika, two viruses never before seen in many of these areas. Furthermore, the latest socioeconomic and political instability in Venezuela and the mass migration of its population into the neighboring countries has given rise to concerns of infectious disease spillover and escalation of arboviral spread in the region. RESULTS: We performed phylogeographic analyses of dengue (DENV) and chikungunya (CHIKV) virus genomes sampled from a surveillance site in Ecuador in 2014-2015, along with genomes from the surrounding countries. Our results revealed at least two introductions of DENV, in 2011 and late 2013, that initially originated from Venezuela and/or Colombia. The introductions were subsequent to increases in the influx of Venezuelan and Colombian citizens into Ecuador, which in 2013 were 343% and 214% higher than in 2009, respectively. However, we show that Venezuela has historically been an important source of DENV dispersal in this region, even before the massive exodus of its population, suggesting already established paths of viral distribution. Like DENV, CHIKV was introduced into Ecuador at multiple time points in 2013-2014, but unlike DENV, these introductions were associated with the Caribbean. Our findings indicated no direct CHIKV connection between Ecuador, Colombia, and Venezuela as of 2015, suggesting that CHIKV was, at this point, not following the paths of DENV spread. CONCLUSION: Our results reveal that Ecuador is vulnerable to arbovirus import from many geographic locations, emphasizing the need of continued surveillance and more diversified prevention strategies. Importantly, increase in human movement along established paths of viral dissemination, combined with regional outbreaks and epidemics, may facilitate viral spread and lead to novel virus introductions. Thus, strengthening infectious disease surveillance and control along migration routes and improving access to healthcare for the vulnerable populations is of utmost importance.

Differential effects of variations in human P450 oxidoreductase on the aromatase activity of CYP19A1 polymorphisms R264C and R264H.

Aromatase (CYP19A1) converts androgens into estrogens and is required for female sexual development and growth and development in both sexes. CYP19A1 is a member of cytochrome P450 family of heme-thiolate monooxygenases located in the endoplasmic reticulum and depends on reducing equivalents from the reduced nicotinamide adenine dinucleotide phosphate via the cytochrome P450 oxidoreductase coded by POR. Both the CYP19A1 and POR genes are highly polymorphic, and mutations in both these genes are linked to disorders of steroid biosynthesis. We have previously shown that R264C and R264H mutations in CYP19A1, as well as mutations in POR, reduce CYP19A1 activity. The R264C is a common polymorphic variant of CYP19A1, with high frequency in Asian and African populations. Polymorphic alleles of POR are found in all populations studied so far and, therefore, may influence activities of CYP19A1 allelic variants. So far, the effects of variations in POR on enzymatic activities of allelic variants of CYP19A1 or any other steroid metabolizing cytochrome P450 proteins have not been studied. Here we are reporting the effects of three POR variants on the aromatase activities of two CYP19A1 variants, R264C, and R264H. We used bacterially expressed and purified preparations of WT and variant forms of CYP19A1 and POR and constructed liposomes with embedded CYP19A1 and POR proteins and assayed the CYP19A1 activities using radiolabeled androstenedione as a substrate. With the WT-POR as a redox partner, the R264C-CYP19A1 showed only 15% of aromatase activity, but the R264H had 87% of aromatase activity compared to WT-CYP19A1. With P284L-POR as a redox partner, R264C-CYP19A1 lost all activity but retained 6.7% of activity when P284T-POR was used as a redox partner. The R264H-CYP19A1 showed low activities with both the POR-P284 L as well as the POR-P284 T. When the POR-Y607C was used as a redox partner, the R264C-CYP19A1 retained approximately 5% of CYP19A1 activity. Remarkably, The R264H-CYP19A1 had more than three-fold higher activity compared to WT-CYP19A1 when the POR-Y607C was used as the redox partner, pointing toward a beneficial effect. The slight increase in activity of R264C-CYP19A1 with the P284T-POR and the three-fold increase in activity of the R264H-CYP19A1 with the Y607C-POR point toward a conformational effect and role of protein-protein interaction governed by the R264C and R264H substitutions in the CYP19A1 as well as P284 L, P284 T and Y607C variants of POR. These studies demonstrate that the allelic variants of P450 when present with a variant form of POR may show different activities, and combined effects of variations in the P450 enzymes as well as in the POR should be considered when genetic data are available. Recent trends in the whole-exome and whole-genome sequencing as diagnostic tools will permit combined evaluation of variations in multiple genes that are interdependent and may guide treatment options by adjusting therapeutic interventions based on laboratory analysis.

Evaluation of 17ß-hydroxysteroid dehydrogenase activity using androgen receptor-mediated transactivation.

17ß-Hydroxysteroid dehydrogenases (17ß-HSDs) catalyze the reduction of 17-ketosteroids and the oxidation of 17ß-hydroxysteroids to regulate the production of androgens and estrogens. Among them, 17ß-HSD type 3 (HSD17B3) is expressed almost exclusively in testicular Leydig cells and contributes to development of male sexual characteristics by converting androstenedione (A4) to testosterone (T). Mutations of HSD17B3 genes cause a 46,XY disorder of sexual development (46,XY DSD) as a result of low T production. Therefore, the evaluation of 17ß-HSD3 enzymatic activity is important for understanding and diagnosing this disorder. We adapted a method that easily evaluates enzymatic activity of 17ß-HSD3 by quantifying the conversion from A4 to T using androgen receptor (AR)-mediated transactivation. HEK293 cells were transduced to express human HSD17B3, and incubated medium containing A4. Depending on the incubation time with HSD17B3-expressing cells, the culture media progressively increased luciferase activities in CV-1 cells, transfected with the AR expression vector and androgen-responsive reporter. Culture medium from HSD17B1 and HSD17B5-expressing cells also increased the luciferase activities. This system is also applicable to detect the conversion of 11-ketoandrostenedione to 11-ketotestosterone by HSD17B3. Establishment of HEK293 cells expressing various missense mutations in the HSD17B3 gene associated with 46,XY DSD revealed that this system is effective to evaluate the enzymatic activities of mutant proteins.

Autism-associated mutations in the CaVß2 calcium-channel subunit increase Ba2+-currents and lead to differential modulation by the RGK-protein Gem.

Voltage-gated calcium-channels (VGCCs) are heteromers consisting of several subunits. Mutations in the genes coding for VGCC subunits have been reported to be associated with autism spectrum disorder (ASD). In a previous study, we identified electrophysiologically relevant missense mutations of CaVß2 subunits of VGCCs. From this, we derived the hypothesis that several CaVß2-mutations associated with ASD show common features sensitizing LTCCs and/or enhancing currents. Using a CaVß2d backbone, we performed extensive whole-cell and single-channel patch-clamp analyses of Ba2+ currents carried by Cav1.2 pore subunits co-transfected with the previously described CaVß2 mutations (G167S, S197F) as well as a recently identified point mutation (V2D). Furthermore, the interaction of the mutated CaVß2d subunits with the RGK protein Gem was analyzed by co-immunoprecipitation assays and electrophysiological studies. Patch-clamp analyses revealed that all mutations increase Ba2+ currents, e.g. by decreasing inactivation or increasing fraction of active sweeps. All CaVß2 mutations interact with Gem, but differ in the extent and characteristics of modulation by this RGK protein (e.g. decrease of fraction of active sweeps: CaVß2d_G167S > CaVß2d_V2D > CaVß2d_S197F). In conclusion, patch-clamp recordings of ASD-associated CaVß2d mutations revealed differential modulation of Ba2+ currents carried by CaV1.2 suggesting kind of an "electrophysiological fingerprint" each. The increase in current finally observed with all CaVß2d mutations analyzed might contribute to the complex pathophysiology of ASD and by this indicate a possible underlying molecular mechanism.

A mutation in MAP2 is associated with prenatal hair follicle density.

Hairlessness is usually a rare trait in pigs; however, in this study, we found hairless (HR) pigs at a relatively high frequency in 1 pig herd. We observed that, the lower hair shaft density of HR pigs could be mainly attributed to the lower hair follicle density, and during the embryonic period, d 39-45 were a critical stage for the formation of the hair follicle. In this regard, d 41 during gestation was a particularly important point. Hair follicle morphogenesis occurring at an early stage of embryo development is similar to humans and mice. Further analyses of association studies based on single-nucleotide polymorphism chip as well as sequence data, mRNA sequencing, immunohistochemistry, and comparative genomics demonstrated that microtubule-associated protein 2 (MAP2) is a key gene responsible for hair follicle density and 1 missense mutation of A-to-G at rs328005415 in MAP2, causing a valine-to-methionine substitution leads to the HR phenotype. Considering the high homology between pigs and humans, our research has some significance for the study of the mechanisms of skin development, hair morphogenesis, and hair loss in humans by showing that the pig may be a more appropriate model in which to study these processes.-Jiang, Y., Jiang, Y., Zhang, H., Mei, M., Song, H., Ma, X., Jiang, L., Yu, Z., Zhang, Q., Ding, X. A mutation in MAP2 is associated with prenatal hair follicle density.

Apnea Associated with Brainstem Seizures in Cacna1aS218L Mice Is Caused by Medullary Spreading Depolarization.

Seizure-related apnea is common and can be lethal. Its mechanisms however remain unclear and preventive strategies are lacking. We postulate that brainstem spreading depolarization (SD), previously associated with lethal seizures in animal models, initiates apnea upon invasion of brainstem respiratory centers. To study this, we assessed effects of brainstem seizures on brainstem function and respiration in male and female mice carrying a homozygous S218L missense mutation that leads to gain-of-function of voltage-gated CaV2.1 Ca2+ channels and high risk for fatal seizures. Recordings of brainstem DC potential and neuronal activity, cardiorespiratory activity and local tissue oxygen were performed in freely behaving animals. Brainstem SD occurred during all spontaneous fatal seizures and, unexpectedly, during a subset of nonfatal seizures. Seizure-related SDs in the ventrolateral medulla correlated with respiratory suppression. Seizures induced by stimulation of the inferior colliculus could evoke SD that spread in a rostrocaudal direction, preceding local tissue hypoxia and apnea, indicating that invasion of SD into medullary respiratory centers initiated apnea and hypoxia rather than vice versa Fatal outcome was prevented by timely resuscitation. Moreover, NMDA receptor antagonists MK-801 and memantine prevented seizure-related SD and apnea, which supports brainstem SD as a prerequisite for brainstem seizure-related apnea in this animal model and has translational value for developing strategies that prevent fatal ictal apnea.SIGNIFICANCE STATEMENT Apnea during and following seizures is common, but also likely implicated in sudden unexpected death in epilepsy (SUDEP). This underlines the need to understand mechanisms for potentially lethal seizure-related apnea. In the present work we show, in freely behaving SUDEP-prone transgenic mice, that apnea is induced when spontaneous brainstem seizure-related spreading depolarization (SD) reaches respiratory nuclei in the ventrolateral medulla. We show that brainstem seizure-related medullary SD is followed by local hypoxia and recovers during nonfatal seizures, but not during fatal events. NMDA receptor antagonists prevented medullary SD and apnea, which may be of translational value.

Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance.

Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.

The HCM-causing Y235S cMyBPC mutation accelerates contractile function by altering C1 domain structure.

Mutations in cardiac myosin binding protein C (cMyBPC) are a major cause of hypertrophic cardiomyopathy (HCM). In particular, a single amino acid substitution of tyrosine to serine at residue 237 in humans (residue 235 in mice) has been linked to HCM with strong disease association. Although cMyBPC truncations, deletions and insertions, and frame shift mutations have been studied, relatively little is known about the functional consequences of missense mutations in cMyBPC. In this study, we characterized the functional and structural effects of the HCM-causing Y235S mutation by performing mechanical experiments and molecular dynamics simulations (MDS). cMyBPC null mouse myocardium was virally transfected with wild-type (WT) or Y235S cMyBPC (KOY235S). We found that Y235S cMyBPC was properly expressed and incorporated into the cardiac sarcomere, suggesting that the mechanism of disease of the Y235S mutation is not haploinsufficiency or poison peptides. Mechanical experiments in detergent-skinned myocardium isolated from KOY235S hearts revealed hypercontractile behavior compared to KOWT hearts, evidenced by accelerated cross-bridge kinetics and increased Ca2+ sensitivity of force generation. In addition, MDS revealed that the Y235S mutation causes alterations in important intramolecular interactions, surface conformations, and electrostatic potential of the C1 domain of cMyBPC. Our combined in vitro and in silico data suggest that the Y235S mutation directly disrupts internal and surface properties of the C1 domain of cMyBPC, which potentially alters its ligand-binding interactions. These molecular changes may underlie the mechanism for hypercontractile cross-bridge behavior, which ultimately results in the development of cardiac hypertrophy and in vivo cardiac dysfunction.